LA-UR-06-2202 
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Convection that develops behind the shock front during the first second of a core-collapse su- 
pernova explosion is believed to play a crucial role in the explosion mechanism. We demonstrate 
that the resulting turbulent density fluctuations may be directly observable in the neutrino signal 
starting at t > 3 — 4 s after the onset of the explosion. The effect comes from the modulation 
of the MSW flavor transformations by the turbulent density fluctuations. We derive a simple and 
general criterion for neutrino flavor depolarization in a Kolmogorov-type turbulence and apply it to 
the turbulence seen in modern numerical simulations. The turbulence casts a "shadow" , by making 
other features, such as the shock front, unobservable in the density range covered by the turbulence. 
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I. INTRODUCTION 
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The leading proposal for the mechanism of core- 
collapse supernova explosions is that the shock first stalls 
and then regains its energy through heatingof the ma- 
terial behind it by streaming neutrinos 0,0] Within 
this paradigm a crucial role is played by vigorous con- 
vection behind the shock front [3, 3- This convec- 
tion, clearly seen in modern multidimensional simula- 
tions [1 H, Q, H , brings the energy deposited by neu- 
trinos in dense regions to the region behind the stalled 
shock front (for a recent overview and references see, e.g., 
fiof). Besides the shock revival, the convection leads to 
asymmetric accretion onto the protoneutron star, pro- 
viding a compelling explanation for the high observed 
pulsar velocities 8]. Finally, evidence for the convection 
comes from observations of SN1987A that indicate ex- 
tensive mixing during the early stages of the explosion 

m 

Turbulent convective motions create a fluctuating den- 
sity field in the post-shock region. Importantly, these 
fluctuations remain long after the shock restarts and a 
successful explosion is obtained 0,0- In particular, they 
persist over the duration of the neutrino burst (~ 10 
sec) and can thus modulate the MSW flavor transforma- 
tions of the neutrinos. The goal of the present work is to 
show that the turbulence seen in the simulations indeed 
leaves an imprint on the neutrino signal, possibly starting 
from about 3-4 seconds after the onset of the explosion. 
This imprint may replace, or combine with, signatures 
of other features already discussed in the literature, such 
as the passage of the front |12j and reverse [Tj| shocks 
through the resonance layer (see also ^3), or effects of 
low-density bubbles . 
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A core-collapse supernova emits both neutrinos and 
antineutrinos, of all three flavors, with energies ^10-20 
MeV. The spectra and luminosities for v e , v e , and v x 
(x = fi,T,p,,f) are in general different. On the way out 
of the star, the neutrinos undergo matter-enhanced fla- 
vor transformations at certain characteristic ( "resonant" ) 
densities, p H ~ 2.4 x 10 3 g/cm 3 (15 MeV/E„)(0.5/Y e ) 
(H or "high" resonance) and pi, ~ 7 x 10 1 g/cm 3 
(15 MeV/E v ){Q.h/Y e ) (Lor "low" resonance) [p. These 
transformations permute the neutrinos in various flavors 
and determine the neutrino spectra detected on Earth. 

To estimate the effect of the turbulence on the per- 
mutation efficiency, we need to know the spectrum of 
the density fluctuations in the turbulence and the re- 
sponse of neutrino evolution to stochastic fluctuations of 
different sizes. The two tasks are tightly coupled. The 
fluctuations relevant to neutrino evolution turn out to 
be smaller than the resolution of the present simulations 
and, hence, must be inferred from the large-scale fea- 
tures seen in the simulations by physical arguments. At 
the same time, the spectrum of fluctuations in a physical 
turbulence is very different from the "^-correlated" noise 
for which analytical solutions for neutrino evolution exist 
ElliE El El ES H3 ■ The appropriate treatment will 
be developed here. 



III. TURBULENCE 

As we will see, the scales of interest are A < 1 5/ sin 2#i3 
km. Unless sin 26 < 10~ 3 , these are much smaller 
than the local scale height r at t > 4 sec. Due 
to Rayleigh- Taylor and Kelvin-Hclmholts instabilities, 
which are clearly seen in the simulations, Kolmogorov 
cascade into smaller scales should develop. Density (just 
like velocity) will then exhibit power-law fluctuations on 
all inertial range scales, that is on all scales between r 
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and viscous cut-off l v ((„<A<r), 

C{k) ee / dx{Sn{0)Sn(x))e- tkx = C k a . (1) 



For the Kolmogorov turbulence the Fourier transform of 
the velocity correlator has a = —5/3. We will consider 
the range — 2 < a < —1.5. 

Eq. J3> implies that a density variation between 
two points separated by a distance A scales as ~ 



fn"'i\ v ^Cnk a . For a > — 1 the variation has 
a strong unphysical dependence on the UV cutoff, oc 
lv ( - a+1 ^ 2 ; and diverges as l v — > 0. This is obvious for the 
5-correlated noise, (5n e (x)6n e (y)) = const x 8(x — y), for 
which a — 0: the small-scale variation is infinite. Even 
if one tries to make sense of the (5-correlated noise by 
imposing an ad-hoc UV cutoff, one has no connection 
with a realistic turbulence. On the contrary, a < — 1 is 
physical as everything is related to the fluctuations on 
the largest scales. Flavor transformations in a medium 
with this spectrum of the fluctuations, to the best of our 
knowledge, have not been treated in the literature [32| . 



IV. EFFECTS OF FLUCTUATIONS: 
MODEL 



A TOY 



We consider a 2-flavor system which goes through a 
"noisy" level-crossing, with a Kolmogorov spectrum. We 
take the Hamiltonian 
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-A cos 26* + A(x) 
A sin 29 



A sin 29 
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with A ee Am 2 /AE l/ , A(x 
50, 9 = 0.1, A(x) =x + A noise 
of the profile is modelled as 



GFn e (x) I '\/2) and set A = 
(x) . The noisy component 
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(x) =F } k 13 cos[kx + 4>(k)], 



fc=i 



(3) 



where F is a normalization factor, /? = a/2 = —5/6 and 
4>{k) are random phases. The parameters are chosen such 
that for F = the evolution is adiabatic (the adiabaticity 
parameter 7 = it A 2 sin 2 29/\ A' \ 1), and the relevant 
range k > 2A sin 2(9 = 20 (see later) is covered. 

Let us compute the survival probability P of a neu- 
trino of a given flavor (for definiteness v e ) traveling from 
x start — —100 to Xend = 100 for different values of the 
noise amplitude F. Since the process is stochastic, for 
each F we repeat the calculation with different random 
phases. The results are shown in Fig. 2] Three different 
regimes are clearly seen: (i) For small noise (F < 10~ 2 ) 
its effects are negligible; P{v e — > 24) is dictated entirely 
by the smooth component of the profile, (ii) As the noise 
is increased, a stochastic scatter appears; from F ~ 0.2 
to 1 the average survival probability - shown by the di- 
amonds - grows as a power law with F. (iii) Finally, for 
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FIG. 1: The electron neutrino survival probability as a func- 
tion of the noise amplitude F (see Eq. For each F, 
the calculation was repeated 66 times, with different random 
phases 4>(k). The points indicate the results of individual 
runs, while the diamonds show the average values. 



F > 1 — 2, the average P(v e — > v e ) saturates to 1/2; any 
value of P(v e — > v e ) in the interval [0, 1] is equally likely. 

The saturation is easy to understand: in strongly fluc- 
tuating backgrounds, the final state expectation value 
{vf\a\vf) is equally likely to point in any direction in the 
flavor space. The averaged final state is described by a 
density matrix diag(l/2,l/2): it is completely depolar- 
ized, all information about the initial state is lost. This 
result, well-known in the case of the 5-correlated noise, 
is expected in general. 



V. METHOD OF ANALYSIS 

Fig. suggests that to understand the impact of the 
density noise on neutrinos it is enough to derive the ex- 
pression for P(y e — > v e ) in the intermediate power-law 
regime. If the answer turns out to be > 1/2, the noise 
completely flavor-depolarizes the neutrino state and the 
true answer is P(v e — > v e ) = 1/2. Given the significant 
uncertainties in the explosion model, the problem may 
not warrant anything more elaborate. 



VI. DERIVATION OF THE EFFECT OF 
FLUCTUATIONS 

To begin, consider a smooth density profile uq{x). Ro- 
tate the Schrodinger equation, i<f>' = H(x)(j), with H{x) 
of Eq. ||2J), to the basis ip( x ) — V(x)<f> which diagonalizes 
the instantaneous Hamiltonian. The full Hamiltonian in 
the rotated basis has the form (see, e.g., |18| for details) 



Tjraass 
H 



— A m —id9 m /dx 
id9 m /dx A m 



(4) 
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where, introducing Ao(x) = Gpno{x) / '\/2, 

A m = [{Acos29-A a ) 2 + A 2 sm 2 29}-^ 2 , (5) 
sin20 m = Asin26>/A m . (6) 

Now let us add a fluctuating density component Sn(x), 
so that SA(x) = Gp5n(x)/^/2. We again rotate the evo- 
lution equation to the instantaneous mass basis. Impor- 
tantly, however, we choose the basic defined not by the 
total matter term, Ao(x) + SA(x), but only Aq(x). We 
have 



it/}' = 



H™ ass (x)+5A 



cos 29 n 
sin 29 n 



sin 26 m 
- cos 29„ 



1>. (7) 



Transitions between the basis states are driven by the 
off-diagonal terms id9 m /dx and 5As'm29 m . Assuming 
the evolution without the fluctuations would be adia- 
batic 33], we neglect the first of these terms and con- 
centrate on the second. Moreover, we assume A m (x) ^> 
8A(x) cos2(9 m , i.e. approximate the diagonal splitting by 
2A m {x). Then, in the perturbative limit, the probability 
of a transition is 



P 



dxSA(x) sin 29 m {x)e^ x dx ' 2A ™ 



Using Eq. Q, we get 



G 2 p 



dk 
2^ 



C(k) 



dx sin 29 m (x) 



exp[i J dx' (2A m (x) — k)] 



(9) 



Let us explore the implications of this result. First, 
consider a turbulent density field superimposed on a con- 
stant smooth density, A (x) = const. After the neutrino 
travels a distance L ^> A" 1 , we have 



P' 



G 2 F L sin 2 29 m C (2A r 



(10) 



The contributions come from a narrow interval of k, with 
the width of ~ L^ 1 , centered on kg — 2A m . 

Eq. I|10|l says that: (i) The depolarization is efficient 
near the resonance, where sin20 m is large and A m is 
small, (ii) The effect is proportional to L and may be 
significant even off-resonance for a sufficiently extended 
turbulent region. We will investigate both in turn. 

Consider now a resonance crossing with a linear 
smooth component, A = A' x. The integral over x in 
Eq. © can be approximately evaluated using a station- 
ary phase method. For each k > 2Asin2#i3 there are 
two stationary points, which on physical grounds (e.g., 
integration over energy, variations in the profile) should 
be added incoherently. We obtain 



P 



Gi 



V2n' 



dkC(k)G ( 



2Asin26> 



(11) 



where the spectral weight G(p) is given by 

G{p)^Q(p-l)p- 1 (p 2 -l)- 1 ' 2 . (12) 

The fluctuations that contribute have wavelengths < 
A.-nE v j(Am} sin2(9), or, for the f/-resonance, < 47 km 
x(E v /l5 MeV)(0.3/sin26»i 3 ). 

For power-law C (k) in Eq. JQ) the integral in Eq. Ijll|l 
can be evaluated analytically, 



pv 



er ~ law ~ G ^°° (2Asin26> 13 ) 



x N, 



(13) 



where the order one numerical coefficient N is given by 
N = ^nT (1/2 - a/2) /[2T (1 - a/2)]. 

For the 5-correlated noise (a = 0) Eq. Ijl3|l gives 



P 



5—corr 



irG F C Q Asm29 13 /V2\n' \. 



(14) 



This agrees with the known result, expanded in the 
same perturbative limit. In the physically interesting 
case of the Kolmogorov spectrum (a = —5/3), we get 
N = v^Fr(4/3)/2r(ll/6) s» 0.84 and 



P 



Kolr. 



0.84G f Co(2Asin26>i 3 r 2/3 /\ / 2|7iol- (15) 
(8) 

Eq. I|15|) reproduces the power-law rise seen in Fig. ^ 



VII. APPLICATION TO THE SUPERNOVA 



Let us apply Eq. {EH to the actual supernova case. 
Four seconds into the explosion (0, Fig. 4a), when the 
shock expands to L ~ 5 x 10 4 km, the density in parts 
of the turbulent region drops to that of the H resonance 
|34| . Estimating |tIq| as ~ n Q / L and using n = V2A/G F 
and Sn 2 L - C (2ir/L) a+1 /{2n\a + 1|], we write the depo- 
larization criterion (P > 1/2, with P from Eq. (1151) ) as 



Sn L /n L > f 0- 3 (Q+1)/2 (AL)-( Q + 2 )/ 2 , (16) 



where / = y/n a /[N\a + l|2 a + 3 ] is an order one factor. 
For Kolmogorov's a = —5/3, we get 

6n L /n L > 0.07 6>^ 3 . (17) 

Fig. 4a of ff| shows that inside the turbulent region 
8nL/riL ~ a few. Thus, the depolarization criterion is 
actually satisfied, and by a large margin. Importantly, 
this result is robust to variations in the details of the 
turbulent spectrum. As a is varied from —2 to —1.5, the 
numerical factor in Ijl7|l varies only from 0.04 to 0.25. 



VIII. OFF-RESONANCE DEPOLARIZATION 

By continuity, depolarization must be present even be- 
fore the density reaches its resonant value [35| . Let us 
estimate at what point the depolarization effect becomes 
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significant. One way is to model the post-shock profile as 
a constant density region of extent L ~ a few x 10 9 cm. 
Using Eq. JUJl with sin26> m = Asin2(9/A m , G 2 F C /2 ~ 
Al/k a+1 = Al{L/2n) a+1 and A m ~ (A - A), we get 

pconst ^ L a+2 A 2 ^2 2 13 Al(A o - A) Q ~ 2 /8vr Q+1 . (18) 

Plugging in = 1.5 x 10 7 eV and L ~ 2.5 x 10 14 eV" 1 
we find that P = 0.1 is obtained when the density is 4 
times the resonant value for sin 2 29i% ~ 0.1; 2 times the 
resonant value for sin 2 2#i3 ~ 0.01; and 1.5 times the 
resonant value for sin 2 20 13 ~ 0.001. To interpret these 
numbers, we consult, e.g., Fig. 2 of 13] where contours 
of constant density are shown as a function of time for 
a one-dimensional model of an exploding supernova. We 
see that for 0.01 < sin 2 20 13 < 0.1 significant (P ~ 0.1) 
depolarization occurs already at t = 3 sec. 

To check the robustness of this result to the details 
of the profile, we can also model the smooth compo- 
nent A a (x) in the post-shock region by a parabola, 
Ao(x) = a + bx 2 /2. In the far-off-resonance limit, 
a - A > A sin 20i3, we get from Eq. © 

pparab ^ ^AW% G 2 Co(a _ AcosM -§ ; (19) 
vb 

where N 2 = 2 a - 1 / 2 v ^T(3/2 - a)/T(2 - a). Eq. JTHJ) 
yields an estimate that is quite similar to the one ob- 
tained using Eq. (JTSJl 3(J. 

IX. EXPERIMENTAL IMPLICATIONS 

Let us briefly list the experimental implications of our 
findings. 

(1) When the densities in the turbulence reach a few x 
10 3 g/cm 3 , assuming sin 2 26*i3 > 10~ 3 , the H-resonance 
should become completely flavor-depolarized. This man- 
ifests itself in a characteristic dip of the average energy in 
the detector, or changing rate of the high-energy events. 
This is qualitatively (but not quantitatively) similar to 
what has been discussed in connection with the passage of 
the front shock 13]. Just like in that case, the signature 



appears either in neutrinos or antineutrinos, depending 
on the sign of the neutrino mass hierarchy. An observa- 
tion of the effect would thus determine the sign of the 
mass hierarchy and place a lower bound on $13. 

(2) The onset of the depolarization would be gradual, 
with the off-resonance effects possibly present up to a sec- 
ond earlier. This may be used to distinguish turbulence 
from the shock effects, as the latter are more abrupt. The 
effect may be even more pronounced for the L resonance. 

(3) The depolarization replaces the signatures of the 
front shock and other features in the density range cov- 
ered by the turbulence. Indeed, flavor-depolarized neu- 
trino is described by a density matrix diag(l/2,X/2), 
which commutes with any Hamiltonian and carries no 
information about density features either before or af- 
ter the turbulence. This extends to the Earth effect 
when the L resonance also becomes depolarized, provid- 
ing yet another possible way to distinguish turbulence 
from other effects at future large detectors. The obser- 
vation about the loss of sensitivity to the front shock 
was also made in an interesting recent paper |3flj . Un- 
fortunatel y. l30l considers an ad-hoc (5-correlated noise, 

following Eaia, m m m i3 . 

(4) The density profile varies between models and, 
moreover, even within a given model Q varies signifi- 
cantly depending on the direction to the observer. Con- 
sequently, there is no unique prediction for when exactly 
the depolarization effects appear first. 

A detailed study of experimental signatures in various 
scenarios (different hierarchies, ranges of $13, explosion 
models) will be published elsewhere. 
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